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Theoretical precision in estimates of the hadronic contributions to 

and (Mz)* 

F. Jegerlehner a , 

a DESY, Platancnallec 6, D-15738, Zcuthcn, Germany 

I review recent estimates of the non-perturbative hadronic vacuum polarization contributions. Since these at 
present can only be evaluated in terms of experimental data of limited precision, the related uncertainties pose a 
serious limitation in our ability to make precise predictions. Besides e + e~- annihilation data also r decay spectra 
can help to get better predictions. Here, it is important to account for all possible iso-spin violations in r-decay 
spectra, from which e + e~ cross sections may be obtained by an iso-spin rotation. The observed 10% discrepancy 
in the region above the p may be understood as a so far unaccounted iso-spin breaking effect. 



1. INTRODUCTION 

Non-perturbative hadronic contributions affect 
electroweak precision observables mainly via the 
hadronic excitations in the photon vacuum po- 
larization (charge screening) which leads to the 
energy dependence of the effective fine struc- 
ture "constant" a(E). Of particular interest are 
a(Mz) (precision physics at LEP/SLC) [1] and 
a M = (g — 2) M /2 which has been measured at 
the unbelievable precision of 0.7ppm at BNL [2]. 
Apart from the electroweak effects (leptons etc.) 
which are calculable in perturbation theory, a se- 
rious problem shows up for the strong interac- 
tion effects (hadrons/quarks etc.) for the calcu- 
lation of which perturbation theory fails. Fortu- 
nately, general principles allow us to evaluate the 
problematic contributions via dispersion relations 
from experimental e + e~- annihilation data repre- 
sented usually in terms of the cross section ratio 



R-y(s) 



er(e + e 



hadrons) 



(1) 



The impact is that the errors of the experimental 
cross section data are now a dominating factor for 
the theoretical uncertainties of electroweak Stan- 
dard Model predictions. Therefore an art has de- 
veloped of getting precise results from measure- 
ments of often very limited precision. The situ- 
ation is also a big challenge for precision experi- 
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ments on er(e + e — > hadrons) as currently per- 
formed by KLOE at DA$NE [3] and BABAR at 
PEP [4] . I should mention that the measurements 
of R are a difficult task: besides the needed par- 
ticle identification and the background rejection 
one has to get 



i? exp (s) = 



Ahad (1 + <*> RC) C norm (s) 
■^norm £ 



(2) 



where Ah a d is the number of observed hadronic 
events, iV norm is the number of observed normaliz- 
ing events, e is the efficiency-acceptance product 
of hadronic events while <5rc are radiative cor- 
rections to hadron production. er norm (s) is the 
physical cross section for normalizing events (in- 
cluding all radiative corrections integrated over 
the acceptance used for the luminosity measure- 
ment) and { s ) = 47ra 2 /3s is the normal- 
ization. In particular this shows that a precise 
measurement of R requires precise knowledge of 
the relevant radiative corrections. 

For the normalization mostly the Bhabha pro- 
cess is utilized [or /i/i itself in some cases]. In 
general, it is important to be aware of the fact 
that the effective fine structure constant a(/z) en- 
ters radiative correction calculations with differ- 
ent scales [i in "had" and "norm" and thus must 
be taken into account appropriately. 

Recent advances/issues in the evaluation of the 
hadronic vacuum polarization effects are based on 
the following results: 

• Updated results from the precise measure- 
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ments of the processes e + e~ — > p — > 7r + 7r~, 
e + e~ — > Cij — > 7r + 7r~7r° and e + e _ — > (f> — > K^Kg 
performed by the CMD-2 collaboration have ap- 
peared recently [5]. The update was necessary 
due to an overestimate of the integrated luminos- 
ity in the previous analysis which was published 
in 2002 [6]. A more progressive error estimate 
(improving on radiative corrections, in particu- 
lar) allowed a reduction of the systematic error 
from 1.4% to 0.6 % . Also some other CMD-2 
and SND data at energies E < 1.4 GeV have be- 
come available. 

• In 2001 BES-II published their final i?-data 
which, in the region 2.0 GeV to 5.0 GeV, allowed 
to reduce the previously huge systematic errors 
of about 20% to 7% [7]. 

• After 1997 precise r-spectral functions be- 
came available [8-10] which, to the extent that 
flavor SU (2)f in the light hadron sector is a sym- 
metry, allows to obtain the iso-vector part of the 
e + e~-cross section [11,12]. This possibility has 
first been exploited in the present context in [13]. 

• With increasing precision of the low en- 
ergy data it more and more turned out that we 
are confronted with a serious obstacle to fur- 
ther progress: in the region just above the u- 
resonance, the iso-spin rotated r-data, after be- 
ing corrected for the known iso-spin violating ef- 
fects, do not agree with the e + e~-data at the 10% 
level [14]. Before the origin of this discrepancy is 
found it will be hard to make further progress in 
pinning down theoretical uncertainties. 

• In this context iso-spin breaking effects in the 
relationship between the r— and the e + e _ -data 
have been extensively investigated in [15]. The 
question remains whether all possible iso-spin vi- 
olating effects have been taken into account in 
which case the discrepancy would have to be at- 
tributed to experimental problems. 

• A new bound <Sa M (0.6 - 2.0GeV) < 0.7 x 
10 -10 [16] for the contributions of mr^f, irrff which 
include decay products from 7r°7, 177, fj, ayy 
yields a severe constraint on possible missing con- 
tributions reported elsewhere [17]. 

• New results for hadronic e + e~ cross-sections 
are expected soon from KLOE, BABAR and 
BELLE. These experiments, running at fixed en- 
ergies, are able to perform measurements via the 



radiative return method [18,3,4]. Preliminary re- 
sults presented recently by KLOE seem to agree 
very well with the final CMD-2 e+e"-data. 

• Last but not least an important change in 
the hadronic contribution to was the change 
in sign of the leading hadronic light-by-light con- 
tribution (ir° exchange) [19]. 

• Progress was made also in calculating the ra- 
diative corrections to 7r + 7r~ production in energy 
scans, for inclusive measurements in radiative re- 
turn [20] and in photon tagging [18] relevant at 
the meson (&,B) factories. 

Some of these results have substantially in- 
fluenced the precision of the evaluations of the 
vacuum polarization effects in OiQED(Mz) and 
(g — 2)^ since 1995 [21]. The present status is 
reviewed in the following. 

2. EVALUATION OF a(M z ) 

The photon vacuum polarization IL^, (q 2 ) mod- 
ifies the fine structure constant according to 

a 



a{q 2 ) = 
Aa = 

n'( 9 2 ) = 



1 - Aa 

-Re (n;(<7 2 )-n;(o)), 



(3) 




and makes it running. The shift Aa is large due 
to the large change in scale going from zero mo- 
mentum to the Z-mass scale n = Mz and due to 
the many species of fcrmions contributing. 

The various contributions to the shift in the fine 
structure constant come from the leptons (lep = 
e, fi and r), the 5 light quarks (u, b, s, c, and b 
and the corresponding hadrons = had) and from 
the top quark: 



Aa = C^had + Aatop 



(4) 



Also ly-pairs contribute at q 2 > 2M^ V (see [22, 
23]). The leptonic contributions are calculable in 
perturbation theory where at leading order the 
free lepton loops yield 



Aaiep(s) 



3^ 

£=e,fi,T 
0.03142 for s 



In (s/mf) - - +0 (mf/s) 



Mf 



(5) 
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where /3g = \Jl — AmJ/s. This leading contri- 
bution is affected by small electromagnetic cor- 
rections only in the next to leading order. The 
leptonic contribution is actually known to three 
loops [24,25] at which it takes the value {Mz ~ 
91.19 GeV) 2 

Aaiep(Mf) ~ 314.98 x 10~ 4 . (6) 

In contrast, the corresponding free quark loop 
contribution gets substantially modified by low 
energy strong interaction effects, which cannot 
be obtained by perturbative QCD (pQCD). As 
already mentioned, fortunately, one can evaluate 

(5) 

this hadronic term Aa^ d from hadronic e + e~- 
annihilation data by using a dispersion relation. 
The relevant once subtracted vacuum polariza- 
tion amplitude (3) satisfies a convergent disper- 
sion relation and correspondingly the shift of the 
fine structure constant a is given by 

A(*W d = -SUL 

3tt \J 4ro 2 s'(s' - s) 

+ f«4=M) <7) 

JE^ ut s'(s'-s) J 

where 

R^s) = 12vrImnL d ( S ) (8) 

is given by (1). Accordingly, the one particle ir- 
reducible (lpi) blob 

which is our relevant building block, is given by 
diagrams which cannot be cut into two discon- 
nected parts by cutting a single photon line. At 
low energies it exhibits intermediate states like 

7T°7, p,U), (/),-■■ , 7T7T, 3-7T, 4lT, • • • , TTTTJ , 7T7tZ, 

• • • , ttitH, • • • , KK, ■ ■ ■ (at least one hadron plus 
any strong, electromagnetic or weak interaction 
contribution). The corresponding contributions 
are to be calculated via a dispersion relation from 
the imaginary parts which are given by the pro- 
duction of the corresponding intermediate states 

2 For m t ~ 174.3 GeV we have Aa t o P (M|) ~ 

d7r lb mi 



in e + e _ -annihilation via virtual photons (at en- 
ergies sufficiently below the point where 7 — Z 
interference comes into play). 

A direct evaluation of the i? 7 (s)-data up to 
— E cut — 5 GeV and for the T resonance- 
region between 9.6 and 13 GeV and applying per- 
turbative QCD from 5.0 to 9.6 GeV and for the 
high energy tail above 13 GeV at Mz = 91.19 
GeV yields 3 : 

A «had( M i) = 0.027690 ±0.000353 (9) 
oT l (Mf) = 128.922 ± 0.049 . 

The contributions from different energy ranges 
are shown in Tab. I 4 . The Euclidean method, de- 
scribed in [30,31], allows to replace data for the 
Adler function by pQCD at space-like momenta 
> 2.5 GeV. This yields 

Aa had(M§) = 0.027651 ±0.000173 (10) 
or 1 {M%) = 128.939 ± 0.024 . 

with a substantially reduced error. This estimate 
is on a sound theoretical basis and should not 
be confused with so called "theory driven" esti- 
mates, which utilize pQCD in a much less con- 
trolled manner. In future this approach would 
allow to evaluate Aa^ d (M^) with an accuracy 
SAa = 0.00007 or 0.00005 provided future cross- 
section measurements allow to reduce the errors 
below Sa<l% up to J/ip or T, respectively. This 
assumes that in the meantime pQCD parameters 
will be known with much better precision as well. 
This reduction of the error by a factor about 5 
is needed in order to satisfy future requirements 
for precision physics at a linear collider [32] . Our 
analysis is as close to the experimental results as 
possible by utilizing the trapezoidal rule together 

3 pQCD for calculating R(s), as worked out to high accu- 
racy in Refs. [26]-[28], is used here only where it has been 
checked to work and converge well: in non— resonant re- 
gions at sufficiently high energies and sufficiently far from 
resonances and thresholds. I have further checked that 
results obtained with my own routines agree very well 
with the ones obtained via the recently published program 
rhad-1.00 [29]. 

4 Table 1 also specifies largely details of the error handling. 
The different energy ranges mark typical generations of ex- 
periments within which systematic errors are considered 
to be 100% correlated, while all errors are treated as in- 
dependent for all entries of the table. 
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^SyblJ 






X PT 


(0.28, 0.32) 


0.04 ( 


0.00) ( 


0.00) 


0.03 ( 


0.00) ( 0.00) 


P 


(0.28, 0.81) 


26.16 ( 


0.24) ( 


0.27) 


24.26 ( 


0.23) ( 0.25) 




(0.42, 0.81) 


3.02 ( 


0.04) ( 


0.08) 


2.75 ( 


0.03) ( 0.07) 




(1.00, 1.04) 


4.74 ( 


0.07) ( 


0.11) 


4.07 ( 


0.06) ( 0.09) 


J/i/> 




11.50 ( 


0.56) ( 


0.61) 


4.06 ( 


0.19) ( 0.19) 


T 




1.27 ( 


0.05) ( 


0.07) 


0.07 ( 


0.00) ( 0.00) 


had 


(0.81, 1.40) 


12.92 ( 


0.13) ( 


0.52) 


11.05 ( 


0.11) ( 0.43) 


had 


(1.40, 3.10) 


27.13 ( 


0.11) ( 


0.60) 


15.75 ( 


0.06) ( 0.37) 


had 


(3.10, 3.60) 


5.31 ( 


0.11) ( 


0.10) 


1.90 ( 


0.04) ( 0.04) 


had 


(3.60, 9.46) 


51.49 ( 


0.25) ( 


3.00) 


8.41 ( 


0.04) ( 0.44) 


had 


(9.46,13.00) 


18.59 ( 


0.25) ( 


1.36) 


0.90 ( 


0.01) ( 0.07) 


pert 


(13.0,oo) 


115.59 ( 


0.00) ( 


0.12) 


1.09 ( 


0.00) ( 0.00) 


data 


(0.28,13.00) 


162.14 ( 


0.74) ( 


3.46) 


73.21 ( 


0.33) ( 0.80) 


total 




277.73 ( 


0.74) ( 


3.46) 


74.30 ( 


0.33) ( 0.80) 



Table 1 



Results for Aa[ 1 5 a ) d (M|) data and Aa { ^ d (-s ) data (V^O = 2.5 GeV). 



with PDG rules for taking weighted averages be- 
tween different experiments as described in detail 
in [21]. 

The most important ingredient of our analysis 
are the e + e _ -data which we described in detail 
in [21] (see also [13]) and the new data which have 
become available since then [33] . The distribution 
of hadronic contributions to Aa had in the e + e~- 
data based approach in shown in Fig. 1. 




2. 5. GeV to 

GeV 




-ac to -2.5 GeV 



Figure 1. Comparison of the distribution of con- 
tributions and errors (shaded areas scaled up by 
10) in the standard (left) and the Adler function 
based approach (right), respectively. 



Our results are in good agreement with other 
recent analyzes 



a4 5 2,(m|) 



0.027680 ± 0.000360 [34] 
0.027690 ± 0.000180 [35] 



fll) 



The corresponding values for the effective fine 
structure constant are: 



a _1 (M. 



128.935 ± 0.049 
128.933 ±0.025 



[34] 
[35] 



(12) 



The numerical agreement does not necessar- 
ily mean agreement between the different ap- 
proaches. 



3. EVALUATION OF a, 



(.9 - 2V2 



The anomalous magnetic moment of the muon 
a M provides one of the most precise tests of the 
quantum field theory structure of QED and in- 
directly at a deeper level also of the electroweak 
SM. The precision measurement of a M is a very 
specific test of the magnetic helicity flip transi- 
tion ipL&nv ipji, a dimension 5 operator which 
is forbidden for any species of fermions at the 
tree level of any renormalizable theory. In the 
SM it is thus a finite prediction which can be 
tested unambiguously to the extent that we are 
able to calculate it with the necessary accuracy. 
For the perturbative part of the SM an impressive 
precision has been reached. Excitingly the new 
experimental result from Brookhaven [2] which 
reached a substantial improvement in precision 
shows a 1.9[0.7] u deviation from the theoretical 
prediction: \a cx v - a* llc | = 221(113)[074(104)[ x 
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10 -11 , depending on whether one trusts more 
in an e + e~-data[r-data] based evaluation of the 
hadronic vacuum polarization contribution [14] . 

Again contributions from virtual creation and 
reabsorption of strongly interacting particles can- 
not be computed with the help of pQCD and 
cause serious problems. Fortunately the major 
such contribution again enters via the photon vac- 
uum polarization which can be calculated along 
the lines discussed for the effective charge. The 
contribution is described by the diagram 




hadrons 

H / N fj 

and is represented by the integral 



had 



3tt 



ds 



ds 



4m* 

iffi« CD (a) K(s) 



(13) 



which is similar to the integral (7), however with 
a different kernel K(s) which may conveniently 
be written in terms of the variable 



1-0 



l+0„ 
and is given by 



1 



Am 2 ^/ s 



K(s) 



(2 _ S 3 )+ ( 1 + * 2 )( 1 + *) 2 



ln(l + x) - x + ■ 



g-tf^ (14) 
(1-x) 



The integral (13) is written in terms of the 
rescalcd function 



K(s) 



3s 



K(s) 



which is bounded: it increases monotonically 
from 0.63 at threshold s = 4m^ to 1 at oo. Note 
the extra l/s-enhancemcnt of contributions from 
low energies in a M as compared to Act. The rela- 
tive importance of various regions is illustrated in 



final state 


range (GeV) 


5a M 


(stat) (syst) 


X PT 


(0.28, 0.32) 


2.14 ( 


0.02) 


( 0.03) 


P 


(0.28, 0.81) 


429.02 ( 


4.95) 


( 5.59) 


UJ 


(0.42, 0.81) 


37.99 ( 


0.46) 


( 1.03) 


<t> 


(1.00, 1.04) 


36.07 ( 


0.50) 


; 0.83) 






8.74 ( 


0.41) 


; 0.40) 


T 




0.11 ( 


0.00) 


; o.oi) 


had 


(0.81, 1.40) 


105.17 ( 


1.18) 


( 3.29) 


had 


(1.40, 3.10) 


56.33 ( 


0.21) 


( 1.47) 


had 


(3.10, 3.60) 


4.06 ( 


0.08) 


{ 0.08) 


had 


(3.60, 9.46) 


14.43 ( 


0.07) 


; 0.75) 


had 


(9.46,13.00) 


1.30 ( 


0.02) 


; o.io) 


pQCD 


(13.0,oo) 


1.53 ( 


0.00) 


1 0.00) 


data 


(0.28,13.00) 


693.22 ( 5.15) (6.83) 


total 




694.75 ( 5.15) (6.83) 


Table 2 



(la.ta 



Results for da' 



Tab. 2 and Fig. 2. The update of the results [21] 
including the recent data from CMD-2 and BES 
II yields 

a had = (694.75 ±8.56) x 1(T 10 . 



The most recent BNL (g^ 
gives (world average) 



(15) 



2) measurement [2] 



x 10" 



a^ xp = (11659203 ± 8 
which compares with the theoretical prediction 5 
(11659169.6 ± 9.4) x lO" 10 (SM) . 



,thc 




1.0 GeV 



Figure 2. The distribution of contributions and 
errors (shaded areas scaled up by 10) for aj 1 1 ad . 



5 Recent new results concern the hadronic light-by-light 
contribution [36] and the 0(a 4 ) QED contribution to 
a e [37]. 
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The new analysis [14] is "data-driven" like [21, 
13] and confirms a substantial discrepancy be- 
tween e + e~- and r-data. The r-based result 
agrees with the corresponding result of [13]. The 
status is illustrated in Fig. 3. We refer to Ref. [38] 
for a recent review and possible implications. 



H 



|CERN (79) 

Theory KNY (85) 
j+ E821 (98) 

•j+ E821 (99) 

>j+ E821 (00) 

World 
T E821 (03) 

E821 goal 



EJ 95 (e+e") 

181.3 ± 16.0 
DEHZ03 (e+e") 

180.9 ±08.0 
DEHZ03 (t) 

193.6 ±06.8 
FJ03 (e+e") 

168.3 ±08.6 
HMNT03 (e+e- 

166.9 ±07.0 



H 



incl.) 



M 
M 



(a„-116590l)0; X 10" 11 



1.9 a Theory (e + e~ 
0.7 a Theory (t) 
1.9 a 



Figure 3. Experimental (upper part) and theo- 
retical (lower part) status of a M . 



4. e+e CROSS-SECTIONS VIA 

r DECAY SPECTRAL FUNCTIONS 

A substantial improvement of the evaluation of 
a had wou i(j kg possible, by including the r-data, 
provided one would understand iso— spin violating 
effects sufficiently well [15]. This has been pio- 
neered by Ref. [13]. Here one utilizes the fact 
that the vector-current hadronic r-decay spec- 
tral functions are related to the iso-vector part 
of the e + e~— annihilation cross-section via an iso- 
spin rotation: 



X~v T 



X° 



where X~ and X° are related hadronic states. 
The precise relationship may be derived by com- 
paring diagrams like: 




U U 

which for the e + e~ 



r(°) 



cr (e + e ► TV TV 
and for the r case into 



case translates into 
47ra 2 



x v (s) (16) 



i dr. 
r^ (r ~ 

B(t~ 



7T tv®u t 



6ir\V u d\ 2 S E w 



1 



1 



x V-{s) x 
(17) 



2s 

r i- j 2 



B(t — > V T 7T 7T°) 

where \V ud \ = 0.9752 ± 0.0007 [39] denotes 
the CKM weak mixing matrix element and 

'S'EW(ncw) = 1-0233 ± 0.0006 [SEW(old) = 

1.0194] accounts for electroweak radiative cor- 
rections [40-42,15]. The spectral functions are 
obtained from the corresponding invariant mass 
distributions. The B(«)'s are branching ratios. 
SU(2) symmetry (CVC) would imply 



-00 = «o(*0 



(18) 



The spectral functions Uj(s) are related to the 
pion form factors F*(s) by 



Vi(s) 



12tt 



|F l ( S )| 2 ; (< = 0,-) 



(19) 



where (3f(s) is the pion velocity. The difference 
in phase space of the pion pairs gives rise to the 
relative factor 8 3 _ /P 3 - +■ 

With the precision of the validity of CVC, thus 
the r-data allow us to improve the I — I part 



of the 



cross-section which by itself is not a 



directly measurable quantity. It mainly improves 
the knowledge of the 7r + 7r _ channel (p-resonance 
contribution) which is dominating in ajj ad (72%). 
After taking into account the known iso-spin 
breaking effects [15] the r-data show substantial 
discrepancies in comparison with the e + e~-data 
(about 10% just above the p-resonance) . This 
issue can certainly be settled by the radiative re- 
turn experiments with KLOE [3] at LNF/Frascati 
and with BABAR [4] at SLAC. In fact prelimi- 
nary results from KLOE are close to the CMD- 
2 results. At present one obtains incompatible 
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predictions for a^ ad based on e + e -data or on 
T-data (see Fig. 3). 

5. ISO-SPIN BREAKING 

CORRECTIONS IN r VS. e+e" 




300 400 500 600 700 800 900 1000 1100 
E (MeV) 

Figure 4. The ratio between r-data sets from 
ALEPH, OPAL and CLEO and the I = 1 part 
of the CMD-2 fit of the e+e~-data. The curves 
which should guide the eye are fits of the ratios 
using 8th order Tschebycheff polynomials. 



Before a precise comparison is possible all kind 
of iso-spin breaking effects have to be taken into 
account. As mentioned earlier, this has been in- 
vestigated in [15] for the most relevant tttt chan- 
nel. For the r version of the pion form factor, 
following from (17) and (19), we perform the iso- 
spin breaking corrections 



rm{s) 



1 



ffEW(old) 

Gem(s) Pl-^o Sew( 



(20) 



with Gem(s) (from [15]) accounting for the QED 
corrections of the r-decays. Final state radiation 
(FSR) is modeled by scalar QED. The resulting 
bare form factor \F„ |?J) =1 (s) compares to the bare 
(vacuum polarization and FSR subtracted) form 
factor |i^r|? e +^j:L\(s), which can be obtained from 

the measured pion form factor \F-n-\(e+e-\(s) by 
subtracting the p — u) mixing effects: 



|i^| 2/=1 ( s ) = |^| 2oxp ( s )/| ( 1+ ' " 



(s u - s) 



(21) 



with s w = (M u — fr^) 2 , e determined by a fit 
to the data: e = 0.00172. In Fig. 4 we display 



l^r|^ _1 (s)/|i^|( e +--->(s) which shows large de- 
viations from the CVC line represented by unity. 

These above corrections were applied also 
in [14] and revealed that they were not sufficient 
to remove the unexpectedly large discrepancy (see 
[14] for details). The only large effect I am aware 
of (of order 10%) which is in the game of the com- 
parison is a possible shift of the invariant mass 
of the pion-pairs in the p resonance region. An 
idea one gets if one is looking at the experimen- 
tal p-mass values, shown in the particle data ta- 
bles [39] ("dipole shape"). If the energy calibra- 
tion of the 7T7r-system would be too low in e + e~ 
measurements or to high in r measurements by 
1% one could easily get a 10% decrease or in- 
crease in the tail, respectively. Since the j o ± — p° 
mass difference as well as the difference in the 
widths T^'°(p — > tttt, Trir'-f) are neither experimen- 
tally nor theoretically established, corresponding 
iso-spin violations cannot be corrected for appro- 
priately. Note that the subtraction of the large 
and strongly energy dependent vacuum polariza- 
tion effects necessary for the e + e~-data, which 
seems to worsen the problem, is properly treated 
in the analysis. 

6. EVALUATION OF a M VIA THE 
ADLER FUNCTION 

In Ref. [30] it has been shown how one can ob- 
tain a better control on the validity of pQCD by 
utilizing analyticity and looking at the problem in 
the ^-channel (Euclidean field theory approach). 
It has been found that "data" may be safely re- 
placed by pQCD at > 2.5GeV. An applica- 
tion to the calculation of the running fine struc- 
ture constant has been discussed in [31,33]. Here 
we consider the application to the calculation of 
a]^ ad . Starting point is the basic integral repre- 
sentation 



oo 1 
had OL f ds f 





x 2 (1 - x) 



c) s/rrifi 37r 



R(s).(22) 



If we first integrate over x we find the well known 
standard representation as an integral along the 
cut of the vacuum polarization amplitude in the 
time like region, while an interchange of the order 
of integrations yields an integral over the hadronic 
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shift of the fine structure constant in the space- 
like domain [43]: 

1 

a l ad = f J dx(l-x) Aa had (-Q 2 {x)) (23) 
o 

where Q 2 (x) = j^m^ is the space-like square 
momentum-transfer or 

In this approach we (i) calculate the Adler func- 
tion from the e + e~-data and pQCD for the tail 
above 13 GeV, (ii) calculate the shift Aa had in 
the Euclidean region with or without an addi- 
tional cut in the i-channel at 2.5 GeV and (iii) 
calculate a£ ad via (23). 

Alternatively, by performing a partial integra- 
tion in (23) one finds 

l 

< ad = J dxx(2-x) (D(Q 2 (x))/Q 2 (x))(2A) 
o 

by means of which the number of integrations 
may be reduced by one. The evaluation in both 
forms provides a good stability test of the numer- 
ical integrations involved. 

Utilizing the most recent e + e~-data we ob- 
tain a result which agrees with the values ob- 
tained by the direct evaluation also in the er- 
ror. Not too surprisingly, as is well known, the 
contribution to a^ ad is dominated by the low en- 
ergy e+e~-data below 1 GeV; here the replace- 
ment of data by pQCD does not reduce the uncer- 
tainty. The reason is hat the pQCD contribution 
replacing the Euclidean Adler function at \J~-i > 
2.5GeV shows a substantial uncertainty due to 
the uncertainty of the charm mass m c (m c ) = 
1.15. ..1.35 GeV. The uncertainty in the strong 
coupling constant a s {M\ ) = 0.120±0.003 is small 
and is not the dominating effect. In contrast 
to other authors which use pQCD for estimat- 
ing R(s) directly, we do not obtain a reduction of 
the error. Of course our cut at 2.5 GeV, which 



we think is all we can justify, is more conserva- 
tive than the 1.8 GeV in the time-like region an- 
ticipated elsewhere. Thus the best value we can 
obtain from presently available e + e~-data alone 
is the result (15). The Euclidean method of cal- 
culating a{j ad will only be useful once the QCD 
parameters will be determined much more accu- 
rately. For recent progress in this direction we 
refer to [44-47] . 

7. SUMMARY AND OUTLOOK 

Evaluations of the hadronic vacuum polariza- 
tion effects based on e+e~-data agree reason- 
ably well between different groups, especially for 
the "conservative" analyzes which rely directly on 
the experimental data. The accuracy typically is 
1.3%. Future precision physics at a high lumi- 
nosity linear collider would require an improve- 
ment of the precision in a(Mz) by a factor of 
about 5 at least. The "theory-driven" analyzes 
which replace data by perturbativcly calculated 
i?-ratios obtain much smaller errors (about a fac- 
tor 2) in general. Obviously much more reliable is 
the Adler function "monitored" evaluation, which 
utilizes pQCD only for the smooth Adler function 
which is much better accessible to perturbation 
theory at sufficiently large space-like momenta 
(\Q\ > 2.5GeV). The reason why this approach is 
not so popular is the fact that the method is much 
more elaborate. Necessary improvements are pos- 
sible only by measurements at the 1% level of the 
hadronic cross-sections up to J/tp or better up to 
the T. Plans for future experiments in this direc- 
tion exist at many places (Novosibirsk, Frascati, 
SLAC, Beijing, Cornell and KEK). The r-data 
are potentially important for improving the a^ ad 
evaluation. However, they can only be utilized 
after appropriate iso-spin corrections. It is likely 
that the observed r vs e + e~ disagreement is a so 
far unaccounted iso-spin breaking effect, due to 
the difference in the physical mass and width of 
the charged observed in r-decays and of the 
neutral p° seen in e + e~ annihilation. 
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